Energy dissipation in DC-field driven electron lattice coupled to fermion baths 
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Electron transport in electric-field-driven tight-binding lattice coupled to fermion baths is com- 
prehensively studied. We reformulate the problem by using the scattering state method within the 
Coulomb gauge. Calculations show that the formulation justifies direct access to the steady-state 
bypassing the time-transient calculations, which then makes the steady-state methods developed for 
quantum dot theories applicable to lattice models. We show that the effective temperature of the 
hot-electron induced by a DC electric field behaves as Tch — C"f{Q/r) with a numerical constant C, 
tight-binding parameter 7, the Bloch oscillation frequency il and the damping parameter F. In the 
small damping limit F/fi — >■ 0, the steady-state has a singular property with the electron becoming 
extremely hot in an analogy to the short-circuit effect. This leads to the conclusion that the dissi- 
pation mechanism cannot be considered as an implicit process, as treated in equilibrium theories. 
Finally, using the energy flux relation, we derive a steady-state current for interacting models where 
only on-site Green's functions are necessary. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.45. Gm 
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I. INTRODUCTION 

Formulating strong-field transport in electron lattice 
has always been one of the most ch allen ging theoreti- 
cal goals in condensed matter physica^^. This is more 
true with today's advanced nano-lithography techniques 
where we can now realize electron lattice driven far from 
equilibrium. Even though this is an old problem, we are 
still in search of a firm theoretical paradigm to approach 
the problem in general. One of the central puzzles is 
dissipation. In equilibrium, the presence of an open en- 
vironment in contact with a system introduces thermal- 
ization, and once the temperature is defined, we often 
use canonical or grand canonical ensemble, and do not 
consider the coupling to the environment as an explicit 
parameter. We naturally ask whether such simplifying 
ansatz may be possible in nonequilibrium, at least for 
steady-state description. 

In nonequilibrium, we do not know such tremendously 
simplifying paradigms to take the role of the environ- 
ment as an implicit parameter which can be hidden in 
thermalization process. On the contrary, the dissipation 
is considered as an integral part of the nonequilibrium 
process, and we need to include the dissipation mecha- 
nism explicitly for sound theoretical description. In some 
systems, however, the dissipation process can be simpli- 
fied. In quantum dots (QDs) under a finite DC bias, 
electron reservoirs coupled to the QD also act as en- 
ergy source/drain and the energy relaxation is assumed 
to happen far away from QD and inside a battery. The 
electrical leads are then modeled as non-interacting reser- 
voirs, and the electron transport is viewed as a trans- 
mission problem- . By taking the open limit, the excess 
energy can be taken infinitely far away from the QD, 
and the problem supports steady-state^. Various quan- 
tum simulation methods have been proposed to study the 
transient behaviors of interacting modela^"^. In the limit 
where a steady-state exists, the nonequilibrium state can 



also be described within the time-independent statistical 
mechanics frameworlP', from which HershfielcP proposed 
the nonequilibrium density matrix 
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with the reservoir energy tak for the source {a — S) 
and drain (D) with the continuum index k. ijr^^. is the 
creation operator of the full scattering state as the solu- 
tion of the whole system of quantum dot and the leads, 
/ia is the chemical potential of the respective reservoirs. 
The scattering state formulation, although conceptually 
appealing, has initially been applied only to limited mod- 
elJiSldue to the difficulty of finding the scatte ring st ates. 
In recent years, several steady-state methodsiiHi^l have 
been developed and have been extended to general mod- 
els. 

Recently, the attention of the field has turned to lat- 
tice nonequilbrium problems. Even at a very stage of the 
field, there have been important findings in the nonequi- 
lbrium processes, most notably that electrons under a 
DC electric field seem to build up internal energy quite 
quickly, reaching a quite different steady-state from equi- 
librium strongly correlated states^^'. One of the most 
popular technique of solving lattice many-body models 
has been the dynamical mean-field theory (DMFT}iSl. 
Its success has been well documented in the description 
of the Mott transition and strong-correlation physics in 
equilibrium. While there have been a fair amount of 
DMFT works to electric-field driven lattice systems, the 
validity of the local approximation is still unconfirmed 
and subject to intense debate. 

There have been numerous attempts to simulate 
nonequilibrium physi cs in lattice models, often through 
isolated HamiltoniansiSinHlII suited for quench dynamics 
of cold atom systems in optical lattice, periodicall y driv en 
system^iSlSO]^ g^^^j some basic dissipation mo del^i^lSilMl 



However, in part due to the numerical difficulties in 
simulating long time-evolution, most of the efforts have 
focused on high-field pheno mena s uch as the dielectric 
breakdown in Mott insulators ^^^l The main emerging 
picture of the calculations is that the external field drives 
the electronic lattice systems into hot temperature, gen- 
erally regardless of the nature of many-body interaction. 
Even though the picture is in agreement between many 
groups, the detailed understanding of the nature of the 
hot electron state and its eventual fate in more realistic 
setting is not available. To gain systematic understand- 
ing of such nonequilibrium state, the dissipation should 
be included in explicit models and their analytic behav- 
ior has to be studied with the damping as a controlled 
parameter. 

One of the goals of this paper is to introduce steady- 
state formulation via scattering state method for lat- 
tice with fcrmion baths and comprehensively analyze the 
model to show that the system possesses many proper- 
ties which are expected of physical systems, for instance, 
consistent picture as the Boltzmann transport theory. In 
the process, an argument will be made that the fermion 
bath model and the steady-state methods are a good min- 
imal system to study nonequilibrium strong correlation 
physics. In the previous papci^^ by one of Authors, the 
fermion bath model under a DC electric field has been 
shown to reproduce the key ingredients as predicted by 
the classical Boltzmann transport theory, and to have a 
stead-state solution. The occupation number as a func- 
tion of mechanical momentum has been shown to have 
the Fermi sea shift by the drift velocity proportional to 
the scattering time and the electric field. Furthermore, 
the DC current has been derived to be consistent with the 
Boltzmann transport result applied to nanostructurea^Zl. 

In this work, we further develop the solution to show 
that the scattering state formulation is applicable, and 
therefore a wide range of new techniques can be de- 
veloped to solve the interacting lattice nonequilibrium 
phenomena. Explicit calculations from temporal gauge 
and the Coulomb gauge with scattering state formula- 
tion show that they are completely consistent with each 
other. The Coulomb gauge enables the time-independent 
formalism, making physical interpretations more trans- 
parent. We calculate explicitly the local distribution 
function, from which we derive that the effective tem- 
perature scales as Tcs = Cjifl/T) with a numerical con- 
stant C, tight-binding parameter 7, the Bloch oscillation 
frequency 51 and the damping parameter F. The effec- 
tive temperature exhibits a singular limit of Teg — >■ 00 
for r/n — >■ 0. This proves that one should not take the 
damping as an implicit process, as treated in equilibrium 
theory. Finally we derive, via the energy flux conserva- 
tion with the Joule heating, a general DC current relation 
as a functional of local Green's functions as an extension 
of the Meir-Wingreen formuleP^ to lattice models, and 
confirm the linear response theory. 

The main text of the paper is organized as follows. In 
Section II, the method introduced in Ref. l26i is further 



developed for Green's functions in the temporal gauge. 
In Section III, we introduce the Coulomb gauge and show 
that the Green's functions in both gauges become iden- 
tical in the long-time limit. In Section III we further 
discuss several important nonequilibrium quantities: the 
local distribution function and the effective temperature 
in A, time-evolution of wave-packet in B, dissipation and 
energy fiux in C, and finally the derivation of the DC 
current in interacting models in D. Appendices provide 
detailed analytic calculations. 



II. TIME-DEPENDENT THEORY WITH 
TEMPORAL GAUGE 

To demonstrate the equivalence of the time-dependent 
temporal gauge to the scattering-state formalism with 
time-independent Coulomb gauge, we start with the one- 
dimensional non-interacting model considered earlier^Sl. 
We study a quadratic model of a one-dimensional s- 
orbital tight-binding model connected to fermionic reser- 
voirs (see Fig. fl]) under a uniform electric field E. The 
effect of the electric field for time i > is absorbed in 
the temporal gauge as the Peierls phase ip{t) = V,t ■ 9{t) 
in the hopping integraP^ 7. Here fl = eEa is the Bloch 
oscillation frequency and 0{t) is the step function. The 
time-dependent Hamiltonian then reads 






^(4df-KH.c.), 



(2) 



with dj as the (spinless) electron operator on the tight- 
binding chain on site £, c]^^ with the reservoir fermion 
states connected to the site I with the continuum index 
a along each reservoir chain of length L. The length 
L is taken to infinity, and the time scale L/vp (with 
Fermi velocity of the chain vp) for the wave to reach the 
end of the reservoir chain is considered larger than any 
other time scales in the problem. As discussed in Ref. |26l 
the Hamiltonian can be diagonalized in each fc-sector as 

^W = Efc^fcW with 



7t, 



.t 



Hk{t) = -2-fC0'S,[k + Lp{t)]d\dk+'^(iaC'^^Cka 

a 
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(3) 



with the fermion operators Fourier transformed to the 
wave-vector basis. Here ed{k) = —2'-fC0s{k) is the tight- 
binding dispersion at zero ii^-field. The reservoir states 
formed by cj.^ acts as an open particle source with its 
chemical potential set at zero energy. The problem can 
be solved with Hq = fffc(O) as the unperturbed Hamil- 
tonian and V{t) — Hk{t) — Hk{0) as the time-dependent 




The lesser Green's function can be simplified in a 
straightforward calculation from Eq. ^ following the 
similar procedures as in Ref. [55] in the long-time limit 
(t,i'>r-i) as 



G,<(t,t') 
where 



ds 



ds'GUt,s)J:<is-s')GUs\n, 

J 

(10) 



FIG. 1; One-dimensional tight-binding lattice of orbital de 
under an electric field E. Each lattice site is connected to an 
identical fermionic bath of {cea} with the continuum index a 
along the reservoir chain direction. In the temporal gauge the 
effect of the electric field is absorbed in the hopping integral 
—7 with the Peierls phase e*^'. 



perturbation, 

V{t) = -27 {cos[A; + ip{t)] - cos{k)} d\dk = v{t)d\dk. 

(4) 
This block-diagonal Hamiltonian is nothing but a reso- 
nant level model coupled to a reservoir, with the level's 
energy oscillating in time^. With the perturbation of 
one-body terms of a finite degrees of freedom, one can 
write the Dyson's equat ion for the retarded and lesser 
Green's functions aa^SSni 



k — ^kfl^ ^kfl^^k 



G< = [/ + GI.V]G<„[/ + VG^], 



(5) 
(6) 



where G^ and G^ are the lesser and retarded Green's 
function matrices, respectively. G^g and G^ q are for 
the non-interacting limit. The multiplication of Green's 
function matrices denotes time integration. 

Following Ref. [551 the retarded Green's function is 

Gl{t,t') = -i0{t - t')e-^\*-''\+'^'^' it' '=°^ik+^'')d'> ^ (7) 

with the damping parameter T — (g^ / L) ^^ Sie^) for 
reservoirs of flat density of states of infinite band- 
width. The local retarded Green's function G\^^{t,t') = 
(27r)-i /_^^ Gl{t, t')dk becomes 

G\^At,t') = -^0(t-Oe-^l*"*'IJo (^^sin^^i^) , 

(8) 
with the zero-th Bessel function Jq{x). Here, the gauge- 
invariant local function becomes a function of only the 
relative time, G'{^^{t,t') = G\^^{t — t'). Fourier transfor- 
mation with respect to the relative time gives 



00 J (22L\'2 



m— — 00 



mn + ir' 



(9) 



by using the Bessel function relatiorplj Jo(2zsin^) 



IT 



E<fs)= / —e-^^'duj, (II) 



with the self-energy Sr from the damping taken as the 
perturbation. Although the above equation has been de- 
rived with the time-dependent Peierls term Eq. Q as the 
perturbation, the same result can be obtained when the 
damping is considered as perturbation in the limit that 
transient terms die out. The local lesser Green's function 
can be computed as G<^{t,t') = {2tt)~^ J^^G<{t,t')dk, 
which again renders the Green's function only dependent 
on the relative time. After changing the dummy variables 
s — t ^ s and s' — /;'—)• s', we have 



GfoM 



ir 

■K 



dt 



duj' 



ds 



ds' 
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i{uj—uj')t — iLj'(s — s')-\-r(s-\-s') J 
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X sm — sm 

2 2 



Again by utilizing the Bessel function relatiorPH 
Jo(Va2 + 62-2a6cosa) = J2^ J™(a) J™(5)e*"" 



G<M = 2zrE/(^ + mf^) / 



^-i{iJ+^)s+rs 



X (7 IT 



47 . ns 



with f{x) — 9{—x), the Fermi-Dirac function at zero 
temperature. From the identitjEH, Jm(2zsin|) — 



-iml7T—a 



)/2 



J2iJi{^)Jra+e{z)t 



G<,(c^) = 2zrV/(c^ + mf}) 



E 



Jei^)Ji-m{^) 



m + iV 



Although the two Green's functions, Eqs. ^ and 
have been reduced to familiar forms resembling spectral 
representation, a clear relation between them is not avail- 
able yet. In the following section, we discuss the scat- 
tering state formalism and find relations connecting the 
retarded and lesser Green's functions. 



III. SCATTERING STATE FORMALISM 

We have learned from Ref. [26| and the above calcu- 
lations that the fermion bath model has a well-defined 
time-independent limit for gauge- invariant quantities 
such as local Green's function. This observation and the 
presence of infinite reservoir states prompt us to consider 
scattering state formalism^^l As depicted in Figs. p[a), 
for each site on the main chain there are infinite degrees 
of freedom coupled from each reservoirs. Therefore, we 
can rewrite the quadratic Hamiltonian in terms of the 
scattering states originating from the fermion reservoir 
chains^!. To adopt the time-independent scattering the- 
ory, we use the Coulomb gauge as shown in Fig.KTa) with 
the Hamiltonian, 



H, 



Coul 



= -lY.^4+ide 



E.c.)-^mdldi 



+ y](^" -i^)cla^eo 



ia 



(13) 



^(cJ„d,-|-H.c.), 



where the static Coulomb potential is applied as a po- 
tential slope to the chain and each reservoirs with their 
chemical potential are raised together with the corre- 
sponding tight-binding sites. 

Since the Hamiltonian is quadratic and the chain is 
coupled to open systems, we can decompose the Hamilto- 
nian in terms of the scattering state operators ipg^ orig- 
inating from the asymptotic state c^^ as given by the 
Lippmann-Schwinger equatiorP^l^ 

1 









-la 



[Hg 



J 



n^ — in — C + ir] 

9 1 

^ ea ~ i^ ~ C + it] 



(14) 



with the Liouvillian operator £A = [-ffcoui,^] for any 
operator and Hg = -(ff/V^) ^^^(cL'^f + H-c). The 
scattering state embodies the openness of the fermion 
reservoirs. As implemented by the infinitesimal imagi- 
nary poles given by irj, once the electrons scatter into 
a reservoir they never come back to the tight-binding 
chain. With the scattering state basis, the Hamiltonian 
is rewritten as 



Hcnu] — 



(15) 



with the Fermi statistics applied separately within each 
^-sector 



(V'LV'to) = /(ec)- 



(16) 



Since the Hamiltonian is quadratic, the scattering state 
operator ipj^ is a linear combination of the original basis 
of dl, and cl, , , and 

t t a' ' 




<AA/# 



= 3 



(&) 



FIG. 2: (a) In the Coulomb gauge for the scattering state for- 
mulation, each tight-binding sites and the associated fermion 
baths are on a potential slope with the potential drop between 
neighboring sites as Q.. (b) As electron moves down the poten- 
tial slope (from f = to ^ = 3 in the figure) , it leaves a trail of 
electron-hole pairs in the fermion reservoirs through particle- 
exchange. Since these electron-hole pairs travel indefinitely 
along the infinite bath chains, the reservoirs act like energy 
drains mimicking inelastic processes. There is an energy flux 
to the reservoirs, but no net particle flux. 



with the expansion coefficient for a fermion annihilat 
operator a given as 



;ion 



Ce.a{a) = {a,i/'lo 



J = 



m-c- 



irj 



(18) 
It is important to realize that, for a quadratic Hamilto- 
nian, the anti-commutation of operators in Cia{a) is just 
a c-number. In fact, this c-number is nothing but the 
retarded Green's function G^i^{ea — t^) between d^i and 
d\. From this argument, the retarded Green's function 
in any time-independent quadratic Hamiltonian is inde- 
pendent of statistics. Note that the scattering state V'Iq, 
originating from the ^-th reservoir has admixture from 
any reservoir states c\,^,. Here, we have put an over-line 
to denote the Green's function in the Coulomb gauge. 
We then have 



4a 



9 



L 



G^,^(e„ - en)di + 



(19) 



This equation can be inverted to express d£ in terms of 
the scattering state basis ipia' as 



de = J2^de'^')^i'- 



(20) 



Similarly as above, Ci^l'a') = {V']/„/,rf^} since the scat- 
tering state operators also satisfy the anti-commutation 



relatioi^{i/']^^,'(/'£'a'} = ^w^aa', and 



(21) 



We will discuss below how G^^, (oj) is explicitly calculated. 
On-site retarded Green's function at the central site ^ = 
can be expressed in terms of the scattering state basis as 



E 



G^,(e^~m)[G^,,{e^,~(.'n)Y 



^o"^""^ ~ L ^ u-e^+m + i7^ 



x({V'fa,'0l/„-}) 



^:^ 



|Go,(e„-m)|^ 



L ■^— ' uj — ec^ -\' to. -\- IT] 

la. 



(22) 



Here the retarded Green's functions appear on both sides 
of the equation, and its self-consistency will be examined 
below. 

The lesser Green's function can be calculated similarly 
as 

G«'(w) = i^—^G'iini{(-a-mn)[G],„,{ec-mQ.)]* 

ma 

= 2zr^GL(u;)[G;„(c.)]*/(c^ + mf7), (23) 

rn 

where each reservoir has its own Fermi energy shifted 
by in and {ipl^ipia) = fi^a)- The above expression is 
quite appealing and physically transparent. With the 
dissipation provided by the particle reservoirs, all elec- 
tron statistics are governed by the Fermi statistics of the 
reservoirs and the effective tunneling between site i and 
the reservoir attached at site m is given by the retarded 
Green's function G^^ . It is noted that we use the infinite- 
band approximation for each fermion reservoirs so that 
any reservoir can provide electrons to any other tight- 
binding lattice sites in principle, and all possible thermal 
factors mix throughout the lattice. 

Now, we turn to calculation of retarded Green's func- 
tions. With the time-independent Hamiltonian, we only 
need to invert the matrix as G££,{uj) ~ [M(w)~"'^]££' with 



[M(w) 



{uj + £il + iT)Su> + 7(5|f_£/|^i, 



(24) 



where the retarded self-energy — zF is attached to each 
site i of the tight-binding lattice with the potential slope. 
Solution to the matrix inversion can be founcPSl as 



Gw (w) 



E 



uj + mil -|- iT 



(25) 



which can be easily verified from M(a;)G (uj) = I. 
Substituting this Green's function into Eq. (23) gives 



the identical result Eq. ( 12 ) as derived from the time- 
dependent temporal gauge. Coming back to the re- 
tarded Green's function, we can easily confirm the iden- 
tity Eq. ( 22 ) from a straightforward calculation after sub- 



stituting Eq. ( 25 ) into Eq. ( 22 ) and by using the contour 



integral and the completeness relation of Bessel func- 
tions. From Eq. (25), it follows that 



Gi+k,i'+k{^) = Gu'i^ + k^), 



from which Eq. ( 23 ) satisfies 



Gi+k,i'+ki^) = Gu'{u} + kfl). 



(26) 



(27) 



In an interacting model, the self-energy is expressed in 
Qr,<,> g^j^j inherits the same property. 



^''i+k,i'+ki^) = S«' (^ + k^)- 



Distribution function 



(28) 



The discussion so far has demonstrated explicitly that 
the dissipative system with fermion baths can be de- 
scribed within the steady-state formalism using the scat- 
tering state basis. One of the central quantities to calcu- 
late is the effective local distribution function, 



/loc(w) = - 



ImGo^oH _Y.AGl,{u)\''j{Lo + m) 



2ImGoo(^) 



He \GQt{^ 



(29) 
where Eq. (22) has been used for ImGoo(w) = 



— r^^ |Gof(w)p. This result takes the same form as the 
ansatz considered in Aron et al^^. 

FIG. [3]ja) shows the numerical evaluation of the local 
distribution function /ioc('^)- For F, il ^ 7, /loc(w) is a 
superposition of small steps coming from the thermal fac- 
tors /(w -I- lU) in Eq. (29) with the envelope following a 
smooth profile similar to the Fermi-Dirac function. Even 
though there is no reason to expect that the nonequilib- 
rium distribution mimics the Fermi-Dirac function, we 
can nevertheless fit the result to the function with an ef- 
fective temperature Teg as shown, despite some deviation 
(see Appendix A for more details). As T^s grows towards 
the finite tight-binding bandwidth 47, the fit becomes in- 
accurate. 

The envelope of the local distribution function plays 
a similar role of the Fermi-Dirac function which dictates 
the abundance of electron-hole pairs available for inter- 
action. In the presence of additional many-body interac- 
tions such as the electron-phonon coupling to local opti- 
cal phonons, the available electron-hole pairs for inelastic 
dissipation are given by the /loc(w) profile and the ef- 
fective temperature Toff in the Fermi-Dirac function will 
play the role of hot electron temperature effectively. 

It is remarkable that T^q seems to approach infinity as 
the damping parameter F becomes smaller. Although it 
may look counter-intuitive at first, this is only the mani- 
festation of the short-circuit behavior where a finite volt- 
age applied across a low resistance conductor induces an 
extremely hot temperature. This is also consistent with 
the numerical calculations with the general conclusion 
that the electron temperature reaches an infinity in closed 
interacting models. 
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FIG. 3: (a) Local distribution function fioc{^) for several 
parameters of the damping F and the Bloch oscillation fre- 
quency n (F, n <C 7). The effective temperature is estimated 
through a fit to the Fermi-Dirac function, (b) The effective 
temperature TaS as a function of f2/F. Up to Q/T ~ 1, Toff 
is well described as a linear function of Q/T. The dashed line 
denotes the analytic expression, Eq. (A4|, derived from the 



low frequency fiodoj) as shown in Appendix A. The blue line 
is from the Kubo formula with T^g = (6/7r^)^''^7(n/F). See 
Appendix B and discussions in section III.D. (c) For larger 
field (f2 = 7 = 1), the steps become more prominent and the 
definition of the effective temperature becomes less robust. 
Nevertheless, the trend in (a-b) continues and at F = 0.2 
even a population inversion happens. 



Numerical fit indicates that T^s is an increasing func- 
tion of n/F for a wide range of (^i, F), although the func- 
tional form eventually deviates from the form as shown 
in FIG. Isjb). For small il/F, the effective temperature 



T, 



off 



Cl 



(30) 



with a dimensionless numerical constant C. This equa- 
tion is one of the key results of this paper. In Appendix 
A, we derive the above linear dependence of (rj/F) and 
approximately estimate that the constant C ^ 2 by an- 
alyzing the w = step in /joc {(^) ■ The effective tempera- 
ture has been previously observed in the momentum dis- 
tribution function and has been speculatecP^ to behave 
as Toff = rt^/r based on the DC conductivity analogy. 
More careful and quantitative analysis now shows that 
the correct dependence is the above relation, Eq. (30 1. 



The Toff relation is even further corroborated with C de- 
rived from the Kubo formula (see Appendix B and dis- 
cussions in section III.D.) as the blue line in FIG. [st 
The Kubo formula result 



T, 



cff 



6 



1/2 



7 



0.77967 (^ 



(31) 



should be exact for the limit fi/F — > 0. The two analyt- 
ical estimates bracket the numerical Toff [see FI G. [3[ b)], 
which shows that the expressions Eqs. (30) and (|31|) 



are 



a reliable approximation for F and 51 up to 7. 

The divergent effective temperature should be taken 
with a caution to interpret in finite bandwidth systems. 
Unlike with the quadratic dispersion relation for con- 
tinuum models22l, the kinetic energy in the single-band 
tight-binding model is always bounded and thus cannot 
give off extremely /loi-electrons to the environment. 

As fl and F become comparable to the bandwidth [see 
FIG. Isjc)], the signature of Bloch oscillation steps be- 
come more obvious and the definition of the effective 
temperature as determined by the shape of the overall 
/loc(w) is not very robust. While the infinite tempera- 
ture in a finite bandwidth system may be questionable, 
the trend observed in FIG.IsFa-b) continues. In the small 
damping limit at F = 0.2, even a population inversion 
happened in the local distribution function. 



B. Time-evolution of wave-packet 

So far, we have seen that the steady-state formalism 
provides a convenient theoretical framework for nonequi- 
librium lattice model of fermion baths. To better under- 
stand the nonequilibrium steady-state, we now look at 
time-dependent quantity indicative of wave-packet drift. 
A steady-state by definition is a time-independent refer- 
ence state where a direct observation of time-evolution of 
a moving particle cannot be made. To confirm that an 
electric charge moves down the potential slope as a func- 
tion of time, we create a hole out of occupied states^ at 
a central position and observe its movement to a different 
position Xh{i) as a function of time in 

= |(4(t)do(0))|2 = |G<(-i)|^ 

(32) 



\{xhi£),t\xhiO),( 




2 

site index (, 



FIG. 4: Observation of wave-packet evolution out of nonequi- 
librium steady-state. Disturbance created out of the steady- 
state travels down the tight-binding ladder as time evolves. 
The amplitude of the wave-packet diminishes due to the de- 
phasing provided by the fermion baths. The curves have been 
off-set for better visibility. 



The lesser Green's function can be easily decomposed 
in terms of scattering states and FIG. |4] shows the wave- 
packet traveling in the direction of the applied field. Am- 
plitude of the observable decays as e~^^* due to the de- 
phasing of electrons from the fermion baths. 

As depicted in Fig. [2jb), electrons travel down the po- 
tential slope in the Coulomb gauge by creating a trail of 
electron-hole (e-h) pairs in the reservoirs. As long as the 
bandwidth of the reservoirs is greater than the poten- 
tial drop D, between neighboring sites (we assumed that 
the bandwidth is infinite in explicit calculations), each 
reservoirs can accommodate an e-h pair with its energy 
matching fi by particle exchange via tunneling. For nar- 
rower reservoir bandwidths, multiple e-h pairs should be 
created to establish a DC current. Since the open reser- 
voirs are of infinite length, the created e-h pairs travel in- 
definitely inside the reservoirs and therefore the fermion 
baths can produce effects similar to the inelastic pro- 
cesses in bosonic baths. We discuss this further in the 
subsequent sections. 



C. Dissipation and energy flux 

We turn to discussions of energy dissipation. The 
Hamiltonian, Eq. (13 1, can be divided into three parts 
as 



_ffs- 



H' 



TB 



H, 



bath 



Hr 



(33) 



with each term representing each line in Eq. ( 13 1, respec 



lively. In the steady-state limit, the energies stored in 
Htb and T^coup are stationary J^(-ffTB) = ^(^mup) = 0, 
as will be demonstrated below. Unlike the case with 
-ffxB and Hconp which are of finite spatial extent, the 
energy flux in i?bath can be non-zero. In the scatter- 
ing theorjBl, the scattering states are formulated in the 



limit that the spatial extent of the scattered wave (Lscatt) 
into the reservoirs is much shorter than the length of the 
reservoir chain {L^catt ^ L). Therefore the scattering 
state represents a solution that the scattered wave con- 
stantly propagating inside the reservoirs without being 
backscattered from the edge of the reservoirs, and quan- 
tities involving the extended states c^^ or cik do not have 
to be stationary in generaP^. 

To be more concrete, we discuss explicit calculations. 
With the fermion baths, a DC electric-field establishes a 
DC current J, as calculated in Ref. [2H1 To investigate 
the effect of the Joule heating, we consider TJtb 



dt 



iH' 



TB/ 



.^Coul, -ffxB, 



^{J) + iigY.{{^ 



(34) 



t+i 



_{)df-n.c. 



-i^i. 



with the current operator J = *7X]f ('^l+i'^^ ~ d^^dg^i) 
within the main chain and q = (1/vL) X]q '^(a- We used 
the steady-state condition for the occupation -^{dldg) — 
0. The first term represents the Joule heating and the sec- 
ond term the energy fiux of electrons from the kinetic en- 
ergy of the main chain into the coupling -ffcoup- Denoting 
the energy flux per site as P, we show that (P) = — 0( J). 
We present the detailed calculations in Appendix C and 
show the equality of the above equation based on the 
scattering-state formalism. 

It can be shown further that J^ (^coup) = 0. As shown 
in appendix B, the energy influx to each of the reservoirs 
is nothing but the Joule heating 



dt 



(/ibath) = 2T ujAioc{^)[fioc{^) - fi^)]duj 

= n{J), (35) 



with the local spectral function defined as 

^ioc(w) = --lmG'oo(w). 

n 



(36) 



The lowercase Hamiltonian h denotes the corresponding 
Hamiltonian per tight-binding site. It nright sound para- 
doxical that the energy in the electronic system is non 
stationary, 

(Here the last equality is from the steady-state current 
taken from Ref.jSH) This is due to the fact that, although 
i^sys governs the electron dynamics, there is another part 
of Hamiltonian which should be included for a closed 
system - the battery connected across the tight-binding 
chain. Since the battery loses its stored charge Q with the 
rate of Q = —{J), the electrostatic energy decrease per 
unit cell of the tight-binding chain becomes J^ (/ibattcry ) = 



— il(J), and the total energy H^ot 
stationary in the steady-state. 



K 



sys 



H, 



battery 



The discussion here again confirms the picture depicted 
in FIG. 2Fb) where the fermion baths act as energy reser- 
voirs while the net electron number flux into the reser- 
voirs is zero. Despite their simplicity, the fermion baths 
through their particle-hole excitations play the role of 
bosonic baths, apart from the boson's explicit dispersion 
relation (with the exception of the Luttinger liquid bath) 
and the physics that might occur from the nonlinear ef- 
fect of the bosonic statistics. 



D. Steady-state current for interacting systems 

From the energy dissipation relations above, we obtain 
the useful formula for the steady-state current. 



merfeld expansiorPSlto derive the linear electrical current 



(J) 



2r 



L0Aioc{^)[floc{^) - f{uj)]duj, 



(38) 



where only on-site Green's functions are needed as in 
Meir-Wingreen formula in quantum dot transporll^Sl. To 
recover the Ohm's law for small field that (J) ex. H. one 
should have that the integral goes as O^ as the leading 
order. This is justified since applying a field of opposite 
direction — f2 should not change the local properties and 
the integral should be of order fl^. This argument can 
be used to analyze the linear response limit, as described 
below. 



The above relation Eq. (38), verified explicitly for the 



non- interacting model in Appendix C, can be extended 
to interacting models. The key identities are steady-state 
conditions 

^ {dl.dea) = and ^ (/ibath) = ^( J) , (39) 

which we expect to hold generally for interacting sys- 
tems as long as the interaction potential does not hold 
infinite amount of energy per site, as in Hubbard model. 
Here we used the spin index a. With on-site interac- 
tion, (hdcr) = ensures zero particle flux into the baths. 
The energy flux equation and the Dyson'e equation hold 
as Eq. (C4| and (C5|, respectively. This immediately 



shows that the current, (38), holds for a wide range of 



interacting fermion bath models. 

Steady-state current derived by Meir-WingreerpSl has 
been widely used in quantum dot calculations. The for- 



mula ( 38 ) can be seen as its extension for lattice mod- 



els with fermion baths. The equation is a functional of 
only local Green's functions. However, it should be made 
clear that while the simplified Meir-Wingreen formula for 
a single quantum model only requires Gqj>,{uj) for the 
quantum dot, both of Gl^^{uj) and G^^{uj) are necessary 
in the lattice models. 



Using the key equation i38h , and Eq. ( 30 ) in the non- 



interacting limit, a linear response limit can be analyzed. 
In the limit of fi/F — >■ 0, the effective temperature is 
expected to be small and, therefore, we can use the Som- 



2T tt"^ 
Jo = e— — r<|fAioc(0) = eC^-^Aioc(O) • O, (40) 



.2-'l' 



and we obtain the linear DC conductivity 

ao = Jo/E - C^^e^a-f^AioM ■ T'K (41) 

Here Aioc(O) is the equilibrium spectral function eval- 
uated at the Fermi energy. With 7ylioc(0) ^ 1, '^ ^ 
m*^ and the scattering time t ^ F^^, we recover 
the Drude conductivity^ni. Comparing this with the lin- 
ear response theory using the Kubo formula, we obtain 
CTo = 2e^a7^/(7rFA/r2 -|- 47^) as detailed in Appendix B. 
Noting that Aioc(O) = (tt^F^ -1-472)-!, we have from 
Eq. (|4l| that C = y/6/^ = 0.7796. This result is shown 



as the blue line in FIG.ptb) in comparison to the numer- 
ically obtained Toff. 

In the interacting limit, the effective temperature ex- 
pression Eq. (30) should be modified. We expect that 
the same form Tcs = C*7*(ri/F*) holds with renormal- 
ized parameters C* , 7* and F*. Then the linear response 
equation becomes ctq = C*^'^e'^aj*'^Aioc{0) ■ (F/F*^). 



IV. CONCLUSION 

In this work, we have reformulated the electron trans- 
port in tight-binding lattice driven by a DC electric 
field using both time-dependent and time-independent 
gauges. The time-independent Coulomb gauge with 
fermion baths leads to the scattering state description for 
steady-state, which makes the calculation and interpreta- 
tion more intuitive. Nonequilibrium quantum statistics 
of quantum dot model, as proposed by Hershfield^, can 
be extended to nonequilibrium lattice as summarized in 
the scattering state expressions. 



Pnoncq 



exp 



-^bath^EaV'LV'fa 



la 



H = J2ie^ - m)iPl^4'ic., 



(42) 
(43) 



with the inverse temperature /3bath of the baths. The 
reservoir scattering states (represented by ipia) are 
shifted by the applied electrostatic potential £Q, and the 
chemical potential is simultaneously shifted with the elec- 
trostatic potential. Therefore, the energy spectra govern- 
ing the dynamics and statistics are different in the above 
expression. The formalism provides a natural framework 
for approximations such as the dynamical mean-field the- 
ory (DMFT). 

It has been shown that the fermion bath model, al- 
though quite rudimentary, produces dissipation mech- 
anism consistent with the Boltzmann transport theory. 
In particular, the steady-state effective temperature in- 
duced by the external field depends quite strongly on the 



electric field and the damping. The effective tempera- 
ture becomes divergent as Tog oc fl/T for small damping 
r versus the Bloch frequency ft = eEa {a is the lat- 
tice constant, and E the electric field). Although this 
might look surprising at first, this phenomenon is sim- 
ply the manifestation of the short-circuit effect. It also 
verifies various numerical calculations with the infinite 
electron temperature resulting in isolated lattice models. 
These findings have fundamental implications in nonequi- 
librium quantum statistics in that dissipation processes 
cannot be implicitly included as thermalization as in 
the Boltzmann factor of equilibrium Gibbsian statistics. 
Through the energy dissipation and the Joule heating 
in the fermion reservoirs, a general DC current relation 
Eq. (38) has been derived for interacting models, as an 



extension of the Meir-Wingreen formula to nonequilib- 
rium lattice systems. The linear response limit has been 
confirmed within this formalism. 

Despite the lack of momentum scattering and explicit 
inelastic processes, the generic features of the fermion 
bath model which are consistent with semi-classical the- 
ory are quite significant. Furthermore, for its simplicity 
the fermion bath model can be used as an ideal build- 
ing block for studying strong correlation effects in lat- 
tice driven out of equilibrium. Particularly, with the 
time-independent Coulomb gauge DMFT can be read- 
ily formulated using the scattering state method ^^"^ * ^^ ! 
It is well-known in equilibrium strong correlation physics 
that electrons undergo collective state when a strong in- 
teraction is present, with some emergent energy scale T* . 
One may speculate that an electric field of order fl ^ T* 
would significantly alter the strongly correlated state. 
However, our study suggests that the dissipation strongly 
interplays with the nonequilibrium condition and non- 
trivial physics may arise even at F < 17 <C T*. Further 
systematic studies are necessary to understand the inter- 
play of nonequilibrium and strong correlation effects. 
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FIG. 5: Fit of /loc(i^) at low frequency u. An analytic ex- 
pression of the effective temperature Tcff is estimated from 
the low frequency part plotted in FIG. 3. The slope of the 
fit is approximated from the step A at oj = as —A/Q,. The 
actual fit gives a somewhat smaller slope than the analytic 
estimate. 



Here we look at the limit of small Q and approximate 
Gqq(O) by the equilibrium Green function 



GooH~'«(w + *r) 



1 - 



4^2 



-,1/2 



(w + ^F)2 



(A2) 



Therefore, the analytic expression for the slope of the fit 
becomes 



A 

n 



F 



ily/A'^^+T^ 



F 



(A3) 



By equating this to the slope of the effective Fermi-Dirac 
function (1 -|- e'^^'^<'")~^ , we obtain 



Tcff 



(A4) 



Note that while the actual numerical fit overestimates 
Toff from the analytic expression due to the high- 
frequency contribution, the overall functional depen- 
dence is quite reliable for F, fi < 7. 



Appendix B: Conductivity from linear response 
theory 



Appendix A: Analytic estimate of effective 
temperature from fiod^) 

Here we analytically justify the relation Tog oc 7(i7/F) 
by considering the low frequency steps as shown in 
FIG. [5] The first step A at w = can be expressed 
as 



A 



r|Goo(o)p 

ImG;;o(0) 



From the Kubo formulcP, the linear conductivity can 
be exactly calculated in the small H. limit. For conve- 
nience, we calculate the current-current correlation func- 
tion in the imaginary-time formalism and then analyti- 
cally continue to the real-frequency in the optical conduc- 
tivitjEIl. For the uniform (q = 0) response function in the 
Matsubara frequency iv, the conductivity is expressed as 



{lm[GU0)-']Y\ (Al) ^M = 



i 1 



y^^vlGk{iuJn)Gk{iuJn 



(Bl) 



k.n 
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Here Vk = 27sin(/c) is the group velocity and the Mat- 
subara Green's function for the electron is given as 

(B2) 
with po{e) =r/Tr-{e^+T'^) ^. Performing the Matsubara 
summation and then the analytic continuation ii^ — > cj + 
ir] for finite uj, we have 



cr(w) 



: / dei / de2Po(ei - <^k)Poi<^2 ~ ^k 
/(ei)-/(e2) 






CJ + ei - 62 



ir] 



(B3) 



Taking its real part and the static limit w — > at zero 
temperature, we obtain the DC linear conductivity 



O-Q 



sin fc 



(B4) 



TT Jo 27r (r2 + 472 cos2 fc)2 

_ 26^072 
7rrVr^"+4^' 

with the restored constants e and a. 

Appendix C: Joule heating and energy flux 



The current expectation value measured at the site (. = 
is expressed as (J) = 27lm(d{(io)- Using the scattering 
state basis, we have 



(dido) 



E 

nkk' 



duj 



E 



^l — k — n'J—k'J — k' — n'J — k' 

{uj + kn- iT){uj + k'il + iT) 

JkJk-ldu! 

{uj + kn- iT)[u! +{k- 1)0 + iV] ' 



by using Eq. (25 1 and the Bessel function identities. We 



have suppressed the argument in the Bessel functions. 
With integration of elementary functions we obtain 



(dido) 



7r(17- 



2iT) E -^M o 

' k ^ 



Jk-l 



2 ^'^ (fc - i)2f72 + r2 



+ 'iXk,k-i 



,(C1) 



with Xnm — TT + tan 



-1 nSl 



-tan" 



-1 mH 



This immediately 



confirms that the current evaluated from the scattering 
state basis matches the result in Ref. [251 

For the operator (P) we need to calculate {{c\ + 
i)do). The c^'-operators are expressed with the scat- 
tering state basis as 



F:t 



1 

7l 



E 



4c. 



L 



^r, 



|g;,,(.<,. - i'n)]'4,. 



e'a' 



(C2) 



A lengthy but straightforward calculation gives 



(C3) 



which confirms the identity (P) = —il.{j). 

For the energy flux into the fermion baths, we examine 



(H, 



bath/ 



— i([-ffsys,-ffbath]) 






(C4) 



First we show that \/L ^ ^^ 



(c!„,d(' 



4cto) = 0. 



From the steady-state condition of -^{dldi) — 0, we 



have {[Hcouhdldi]) = -^{dl^-^^di - djd^+i + 



rf]_irff 



dldf 



- -^J^ai^adi - d\cia) = 0. The first term 
is the total flux into the £-th site due to the current 
along the TB chain, and in the steady-state it is zero. 
Therefore we have zero particle-flux into the reservoir, 
V L^^'Y^^{c'i^di — d\cta) =0. The remaining summation 
\/L-^ J2a''a{clade " ^iQa) IS the energy flux measured 
with respect to the ^-th reservoir chemical potential level 
and it should be independent of i. Setting ^ = 0, we can 
rewrite the expression as follows. 

Consider G<^{t) = i{cljO)do{t)} and G<,{t) = 

i{dQ{—t)coa{0)). For the energy flux per reservoir, we 
can write 



{h 



bath/ 



^Y.'o.j[GU^)-G<, 



[uj) 



duj 
2^' 



From the Dyson's equatiorPS, 

'S{u - e„){Go< (c.) - /(c.)[GSo(w) - GS„(c.)]} 



—2'!Tig 



47r2g 



^{^ - <^a)Aoc{^)[floc{^) ~ /(w)] 



(C5) 



2r / c^Aioc(c^)[/ioc(^) - f{uj)]duj. (C6) 



Then 

(/ibath) 



Taking the time-derivative of the above Green's func- 
tions at i = 0, we have from Eq. (C5), 



= 2r 



E4[Giw-G<,(i)] 



a;Aioc(w)[/ioc(w) - f{oj)]duj = (//bath)- 



This can be equated to 



-w 



^{cla[Hcou\, do] + [Hcou\, dl]coa) 

a 

= -ig-f{cl{di + d-i) - {d\ + dl^)co} 

= -J57((cI+c^)do-4(ci+c_i)), (C7) 
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where the translational invariance of the steady-state has been used. Using Eq. ( C3 ) , Eq. (351 is confirmed 
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